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Abstract 



We use QCD sum rules to compute matrix elements of the AB = 2 operators 
' appearing in the heavy-quark expansion of the width difference of the Bs mass 

eigenstates. Our analysis includes the leading-order operators Q and Qs, as well as 
^ ■ the subleading operators R2 and R3, which appear at next-to-leading order in the 

' l/jTif, expansion. We conclude that the violation of the factorization approximation 

for these matrix elements due to non-perturbative vacuum condensates is as low as 
1-2%. 
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1 Introduction 



The phenomenon of flavour mixing has been intensively investigated over the last decades. 
The standard model of particle physics provides us with a parameterization of flavour 
physics which is compatible with all data taken up to now. However, we are still lacking 
a fundamental theory of flavour, explaining the three-family structure, the masses and 
mixings and CP violation. 

The phenomenology of flavour mixing has a few peculiarities. In the standard model 
the only source of flavour mixing originates from the "mismatch" between the two mass 
matrices for the up and the down quarks, which is encoded in the relative rotation between 
the eigenbases of these matrices given by the CKM matrix. The mass matrices are induced 
by Yukawa couplings to the Higgs particle, which hints at a relation between electroweak 
symmetry breaking and the origin of flavour. 

CP violation in the standard model is related to an irreducible phase in the CKM 
matrix, which can appear for at least three generations [1] . Putting aside the still unsolved 
mystery of strong CP violation [2], this leads to a few interesting conclusions which are 
confirmed by observation. One of these conclusions is the strong suppression of (CP 
violating) electric dipole moments of quarks and leptons, which is compatible with data. 
However, in a generic parameterization of "new physics" contributions it is hard to avoid 
electric dipole moments exceeding the experimental limits by orders of magnitude. 

A further peculiarity of the standard parameterization of flavour physics is the sup- 
pression of "flavour changing neutral currents" (FCNC's) by the GIM mechanism |3], 
which has its root in the unitarity of the CKM matrix. In particular, FCNC processes 
with AB = 2, AS = 2 have been intensively investigated, while AC = 2 processes have 
not yet been observed, in accordance with the very strong GIM suppression predicted by 
the standard model. 

Especially in the systems of neutral B mesons the theoretical description is simplified 
by the fact that the mass difference in these systems is dominated by the short distance 
contribution of the top quark. Furthermore, the width difference, which is expected to be 
sizable in the Bg system, can be computed in the heavy-quark expansion [1]. 

The width difference AP between the Bg mass eigenstates is determined by the off- 
diagonal matrix element of the AB = 2 transition operator T through AP = — 2Pi2 
where 

Pi2 = -^—{BJTlBs) (1) 

and Mb^ is the Bg meson mass. The AB = 2 transitions are initiated by a flavour changing 
neutral current and occur only at the loop level in the standard model. Therefore the 
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transition operator T is a complicated, non-local object. The main problem however 
is the treatment of mesons as bound states of QCD, which involves dynamics in the 
infrared strong coupling regime, where a perturbative treatment is not possible. In the 
heavy-quark expansion the off-diagonal matrix element can be expanded as a series 
in inverse powers of the 6-quark mass as 

mrm = Y.^mo^^='m (2) 

n 

where the Wilson coefficients C„ are calculable in perturbation theory |5j. In this formu- 
lation all the non-perturbative physics is contained in the matrix elements of the local 
AB = 2 operators O^^^^. At leading order in l/nib the transition operator T involves 
two four-quark operators 

Q = ibiSi)v-A(bj^3)v-A (3) 

Qs = (biSi)s-p(bjSj)s-p (4) 

with i a color index. The notation is such that (biSi)v-A = ^i7/i(l — 75)^1 and (biSi)s-p = 
bi{l — 75)sj. At next-to-leading order in l/mb the transition operator involves five new 
(subleading) operators. The complete list of subleading operators and different choices of 
basis can be found in [6l [7] . We shall focus on the operators involving an extra covariant 
derivative acting on the strange-quark field, of which there are four. Neglecting higher- 
order terms in the l/rUb expansion these can further be reduced to the two operators 

R2 = -\{bilD ^D''si)v-A{biSi)v~A (5) 
mi 

R3 = -^{h'D^D^'si)s-p{hsi)s-p (6) 
mi 

with = d^ — igsA^ the covariant derivative. The subleading operators should be under- 
stood in HQET even though they are written formally in terms of full QCD fields. This 
means that the covariant derivative acting on the 6-quark field in ([S])-® can be replaced 
by m^f with v the velocity of the heavy 6-quark, making explicit that the subleading 
operators R2 and -R3 are suppressed only by one power of l/m^. 

The standard parameterization of the matrix elements of these operators is obtained 
through the vacuum saturation approximation [8] with bag parameters Bi controlling the 
accuracy of the factorization, {Bs\Oi\Bs) = Bi{Bs\Oi\Bs)^°''^. For the operators considered 
here, we have (e.g. [6]) (we now use B for the B^ meson and also for the bag parameter 
of the operator Q) 

{B\Q\B) = /|M|2(l + i-)i? (7) 
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(9) 



where A^c = 3 is the number of colors in QCD and Jb is the Bg meson semileptonic decay 
constant. 

The dominant theoretical uncertainties in the prediction of AF = — 2ri2 using the 
heavy-quark expansion are related to the hadronic matrix elements of the local operators 
Oi G {Q,Qs, R2, R3}, or equivalently, the bag parameters Bi. The calculation of the 
bag parameters involves strong interaction dynamics in the infrared region and is thus 
a problem in non-perturbative QCD. The ultimate solution can be provided by their 
direct calculation in lattice QCD. Results for Q and Qs are available, although not yet 
completely reliable [9]. However, a computation for the operators R2 and -R3 is completely 
lacking, and to match the increasing precision of the experimental data it is necessary to 
consider deviations from Bi = 1 even for these subleading operators [7J. 

In this paper we use the technique of QCD sum rules to provide a first estimate of the 
bag parameters for the subleading operators R2 and R^. We focus on the calculation of 
the parameters ABi = Bi — 1, which measure the deviations from the factorization result 
Bi = 1. We limit our analysis to the non-perturbative vacuum condensate contributions 
to these quantities. While more sophisticated treatments with lattice QCD exist for the 
leading-order operators Q and Qs, and with QCD sum rules for Q, we also include these 
operators in our analysis. Studying the full set of operators simultaneously helps clarify 
the general features of sum rules as applied to this class of matrix elements. 

Our main finding is that the non-perturbative contributions to AB are quite small for 
each of the four operators, no larger than 1-2%. We use a simple analytical analysis based 
on the HQET limit within finite energy sum rules (FESR) to give insight into this result. 
To explore corrections to the HQET limit and to provide error estimates we perform a 
more thorough numerical analysis using Borel sum rules. The numerical results suggest 
that corrections to the HQET limit may be large in some cases. 

The paper is organized as follows. In Section [2] we describe the technique of sum rules 
as applied to our case and introduce some necessary notation. In Section [3] we describe 
the calculation of operator-product expansion (OPE) expressions for the Green functions 
used in the analysis. Sections H] and [5] contain our sum-rule analysis and includes full 
QCD and the HQET limit in FESR and Borel form. In Section [6] we give the final results 
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and discuss the assumptions made and uncertainties involved. In Section [7] we give the 
summary of the paper. Some long formulae for the OPE spectral densities are collected 
in the Appendix. 



2 Sum rule calculation of the bag parameters: 
the technique 

In this section we review the sum-rule method for calculating the hadronic matrix elements 
of the AB = 2 operators. The starting point is the three-point correlator 

T{pi,P2) =1^ j d'xdSe'^'--'^'y{Tj{x)0{Q)j{y)). (11) 

The operator O E {Q, Qs, -^2,-^3} is a generic four-quark operator and the interpolating 
current j for the 5-meson can be either an axial-vector (AV) current or pseudoscalar (PS) 
current, defined as 

= s-i^'^^h (AV current) (12) 
js = si'-^^h (PS current). (13) 

The overlap of the interpolating currents with 5-meson states is defined through the 
matrix elements 

(oi^7m75&(o)|5(p)) = if BP,, mi.bmm) = (14) 

where Jb is the semileptonic decay constant of the B meson, Mb is the 5-meson mass, 
rrify is the 6-quark mass, and is the strange-quark mass. For the axial-vector interpo- 
lating current the three-point correlator is a tensor, and we focus on the scalar function 
multiplying the tensor structure PiP2- 

T^''{puP2) = f d'xd'ye'P^'-^-y{Tj^{x)O{0)f,{y)) = p^^T(pi,p2) + • ■ • (15) 



where the ellipsis denote other tensor structures such as PiPi, ^2^2? P1P2 9fj.v It is 
convenient to use the dispersion relation 

dsids2-, 2\ (16) 

(Si -Pl)[s2-Pi) 

and work with the spectral density p(si, 52, g^). Here q = Pi — P2 and = at the 
physical point relevant to the mixing. To derive the sum rules the spectral density is 
evaluated in two ways: 
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1. In a phenomeno logical hadronic picture. In this case the spectral density is modeled 
by a i?-meson pole plus a continuum contribution. This yields 



p\visi,S2)= fl5{s, - Ml)5{s, - Ml){B\0\B) + p^^ 



(17) 



for the axial- vector current, and 



Pps(^1,^2) 




-5{s,- Ml)5{s,~ Ml){B\0\B) + p'^^^' (18) 



(mfc + nis) 



for the pseudoscalar current. 

2. With QCD using the operator-product expansion. The resulting spectral densities 
pP^^ are the sum of a perturbative contribution and a non-perturbative contribution 
involving the vacuum matrix elements of local QCD operators (condensates). 

The idea of QCD sum rules is to use duality between the physical spectrum measured 
in terms of hadrons and the OPE prediction expressed in terms of quarks and gluons 
(the degrees of freedom of the QCD Lagrangian). Duality is implemented by comparing 
integrals of the two spectral densities 



It is common practice to model the continuum contribution to the hadronic spectral den- 
sity with the theoretical expression from the OPE. We choose to match the two expressions 
at the point si = S2 = sq, so that the integration region A in the duality integral is the 
square < Sj < Sq in the (si, S2) plane. One then obtains the sum rules 



Calculating the OPE expressions for the spectral density thus allows for the extraction 
of the hadronic matrix elements {B\Oi\B), or, equivalently, the bag parameters Bi. The 
sum-rule results depend on the parameter sq at which the hadronic continuum is modeled 
by the OPE result; we shall discuss different ways of choosing this parameter later on. 

The sum rules f l20l [2T]) are referred to as "finite energy sum rules" (e.g. [lOj). It is 
expected that results obtained with these basic sum rules give a reasonable approximation 
to a more sophisticated analysis. However, it is also useful to consider a different averaging 
procedure in the duality integrals. The most popular technique is the Borel sum rule 




(19) 




(AV current) 



(PS current). 



(20) 



(21) 
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analysis. In Borel sum rules one works with duality integrals of moments of the spectral 
densities rather than with the spectral densities themselves. In particular, one compares 
the derivatives d^/{dp'^)^ of the spectral densities for large n. In the limit n ^ oo and 
again modeling the hadronic continuum with the OPE prediction one arrives at the Borel 
sum rule 



fl{B\0\B)e~^~^ = f dsids2e~^~^ p^^^{si,S2) (AV current) (22) 

J A 

and analogously for the pseudoscalar case. In the Borel sum rule contributions from 
excited states are exponentially suppressed. Also, studying the stability of the sum rule 
results under variations of the Borel parameters Mi and M2 helps assess their reliability. 

The procedure sketched above can be used to compute the bag parameters directly. 
However, at the level of the OPE, one can identify the contributions to the three-point 
correlator which lead to the value B = 1 only [TTl [12]. Such contributions can be ex- 
pressed as the product of two color-singlet two-point functions, each depending on a single 
momentum. Subtracting this trivial part from the QCD sum rule allows us to focus on the 
piece responsible for deviations from the factorized value. We thus split the three-point 
correlator into two pieces according to 

T{PUP2) = Tfac(pi,P2) + AT(pi,p2), (23) 

where the sum rule obtained from the factorized piece Tfac yields B = 1. This factorized 
part has the explicit form 

TfaciPuP2) = const X n(pi)n(p2) (24) 

with the "const" and the Il{pi) specific to the operator involved. For instance, for the 
operators involving a V-A Dirac structure, one has 

Tt:cipi,P2) = 2 (i + ^) n^(pi)n^(P2) (25) 

with 

p^U^{p) = i J dxe'P^Tj{x)h''{l - 75)s(0)). (26) 

Using this same notation for the factorizable and non-factorizable contributions to the 
spectral densities one finds a sum rule for AB = B — 1 directly. It reads 



flAB{B\0\B)f'''e ^ ^ = J dsids2Ap'l^^{s,,S2)e ^ ^ (27) 
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for the Borel sum rule with an AV interpolating current and analogously for the other 
cases. 

If AB is numerically small compared to the factorized value B = 1 (as one expects 
from the previous analyses [121 lEl [E] and the present study confirms), then this setup 
allows for an essential improvement in precision in comparison with the analysis of the B 
parameter itself. 

2.1 Sum rules in the HQET limit 

The AB = 2 operators are identified by evaluating the transition operator as a series 
in l/m^, according to the heavy-quark expansion. In this treatment, the operators are 
defined in terms of QCD fields and contain implicit mb dependence. For processes contain- 
ing heavy quarks it is advantageous to make this dependence explicit by performing 
calculations in the formal limit — > oo using the framework of HQET. The effective 
theory sets up a systematic expansion in powers of I /mi,, and separates the perturbative 
effects occurring at the scale rrib from those responsible for the hadronic dynamics at the 
scale Aqcd- In addition to our QCD results, we shall consider our results evaluated in 
the HQET hmit. 

To carry out this expansion to a given order in ^^(mfe) and l/ruf,, one must match the 
interpolating currents and the QCD Lagrangian onto their HQET expressions, and evalu- 
ate the three-point correlator in the sum-rule analysis using these effective-theory objects. 
In this paper we shall limit the HQET expansion of a given matrix element to leading 
order in both perturbative and 1 / mi, corrections, ignoring even the effects of leading- log 
resummation. To this level of accuracy the matching onto HQET is trivial, and can be 
obtained directly from the QCD sum-rule expressions by making certain substitutions 
and then expanding in a series in the large 6-quark mass. On the phenomenological side 
of the sum rules, this is done by writing Mb = mi, + A and expanding to leading order in 
A/ mi,. On the OPE side, this is done by writing the spectral variables as Sj = {mb + Ei)"^ 
and expanding to leading order in Ei/mb- 

Applying the HQET expansion to the finite energy sum rules (120112 ip is straightforward, 
and will be discussed in Section 14. 2[ In our numerical analysis in Section [5l we will also 
need the HQET limit of the QCD Borel sum rule fl22l) (and its PS analogue). To obtain 
the HQET expression, we choose the Borel parameters = M| = and define 
W = M^/mi,. Performing the HQET expansion yields 

fl{B\0\B)uQET = 4 / dEidE2e'^''w'''''' p^^l^i, ^2) (AV current) (28) 

J A 
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(a) (b) 

Figure 1: The leading-order perturbative contribution to the three-point correlator (a), 
and a non-factorizable perturbative contribution at next-to- leading-order (b). 

where the HQET limit of the matrix elements are defined by the expansion of the right- 
hand side of d?]). In this case the duality interval A is given by < Ei < Eq. The 
expressions for AB are then derived as before. 



3 The OPE for the three-point correlators 

In this section we describe the calculation of the spectral density functions using the 
OPE. The leading-order results are given by the bare quark loops shown in Figure [U^a). 
The cross denotes the insertion of any one of the four-quark operators Q, Qs, R2, R3, 
and the solid dots can be either axial-vector or pseudoscalar interpolating currents. The 
analysis works very much the same for each of these eight possible cases. Corrections to 
the leading-order result come from two sources: higher-order perturbative corrections and 
non-perturbative corrections in the form of vacuum condensates. Our focus in this paper 
is on the vacuum condensate contributions, which we consider up to dimension six by 
calculating the gluon condensate, the mixed quark-gluon condensate, and the four-quark 
condensate. 

The leading non-perturbative contributions involve the gluon condensate, a dimension- 
four object defined through the vacuum matrix element 

b ^""^ 

{G%G^f^) = ]^2(Ar2 - i^ ^3iiagiyf3 - g/^f^gua) {GG) . (29) 

The non-factorizable corrections proportional to the gluon condensate are obtained by 
calculating the diagram shown in Figure [2|^a) along with the three other permutations 
where the gluons are attached to different loops. Diagrams where the two gluons are 
attached to the same loop are factorizable and hence do not contribute to AB. 

The calculation is most easily performed using the external-field method [15]. The 
advantage of this technique is that the external gluon field can be expressed in terms of 
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(a) (b) (c) 

Figure 2: Non-factorizable contributions involving (a) the {GG) condensate, (b) the {sGs) 
condensate, and (c) the (ssss) condensate. 

the field-strength tensor according to the relation 

A^^{x) = ^x"Gl^ + 0{x'). (30) 

This property allows for a direct extraction of the gluon condensate contributions from 
the diagrams in Fig[2](a), and also simplifies the calculation for the subleading operators 
R2 and R3. Since the operators R2 and are evaluated at the point x = 0, the diagrams 
where a gluon is emitted from the operator itself (the cross in the diagrams) vanish, and 
one need only consider derivative couplings, whose evaluation is essentially the same as 
for the leading-order operators Q and Qs- 

We next consider the dimension-five contributions. These are proportional to the 
mixed quark-gluon condensate, which is defined through the matrix element 

(s^tgsG'^^^t^s,) = ^^^§^{sGs). (31) 

The relevant non-factorizable diagrams are shown in Figure [2](b). As with the gluon 
condensate, the relation (1501) leads to simplifications for the subleading operators i?2,-R3- 
Also in this case one need not consider gluons emitted from the covariant derivative; 
moreover, the external strange-quark fields carry vanishing momentum, so derivatives 
can only act on the strange-quark field contracted inside the loop. 

Finally, we consider the dimension-six contributions involving the four-quark conden- 
sate. The relevant non-factorizable diagrams are shown in Figure [2](c). These vanish for 
the subleading operators R2 and R3, as can be seen by using fl30|) and then noting that 
the derivative terms act on the vacuum fields and thus vanish. For the leading-order op- 
erators Q and Qs the contributions involve matrix elements of the form (srissr2s) where 
the Fj involve both Dirac and color indices. To evaluate these non-factorizable four-quark 
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matrix elements we use the vacuum saturation approximation, by which the full matrix 
element is expressed as 

(sFiSsras) = — ^ (trli tiT^ - trriFa) (ss)^. (32) 

This approximation dates back to the first applications of the sum rule method [16], and 
since then has been checked through numerical analysis in many physical channels. One 
particular study for vector-vector and axial-axial channels established that the factor- 
ization is accurate within 15-20% [17j. Upon using this approximation for the current 
correlator, we find that non-factorizable contributions from the four-quark condensate to 
sum rule for Q and Qs also vanish. Details are given in the appendix. 

We shall limit our OPE analysis to these non-perturbative condensates. To this level 
of accuracy, the OPE result for the spectral density can be written as 

Ap,(si, S2) = App(si, S2){GG) + Apf^s,, S2){sGs) + ... (33) 

for each of the eight cases. Explicit results for the Api can be found in the appendix. 
The result for the operator Q with an AV (PS) interpolating current was first obtained 
in [12] ([H]), while the others are new. The ellipsis refers to the corrections not taken into 
account in our analysis. These include contributions from the dimension six condensate 
(^f'^bcQa^Qb^Qc^\^^ whose numerical value is considered to be small [16]. An attempt to 
take into account condensates of operators of dimension 7 and even 8 was made in ref. [13] 
for the operator Q. We note, however, that the numerical values of these condensates are 
very uncertain and their effects small, and thus exclude them from the analysis. 

More important are higher-order perturbative corrections. The next-to-leading order 
corrections are parametrically on the order of as(mb)/7r ~ 0.06 for as{mi,) = 0.2. Non- 
factorizable perturbative corrections require the evaluation of three-loop diagrams such 
as that shown in Figure [Tj^b). These were calculated in [18] for the leading-order operator 
Q, but are unknown for the other cases. 

As an example and to introduce notation we give here the explicit expression for the 
Qs operator with a pseudoscalar interpolating current: 

Apps(si,S2) = -^{—GG)^— ('^^{6 -3zi-3z2 + Z1Z2) + {piP2fziZ2] 
487r^ 71 S1S2 V 2 / 

((-2 + z,)6{s2 - ml) + (-2 + ^2)^(^1 - ml)) (34) 
167r^ \ / 

Here Zi = ml/si, and 6{si — ml) is the Dirac S function. At the physical point = 
(Pi — P2)^ = the scalar product (P1P2) should be understood as (P1P2) = (^i + ■S2)/2. 
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We also need the HQET expansion of the spectral density, which we obtain by using 
Si = {mh + EiY and expanding to lowest order in Ei/rrib- In this case this limit reads 



4 The bag parameters from finite energy sum rules 

In this section we present the sum-rule results for the ABi using the finite energy sum 
rules ( 1201 [2T|) evaluated at leading order in the HQET approximation. We first give simple 
analytical expressions for the AB^, obtained by relating the sum-rule parameter sq to the 
5- meson decay constant /b, thereby eliminating one parameter. Upon inserting numerical 
values it becomes clear that AB is suppressed by a small scale ratio, independent of the 
particular operator being considered. 

4.1 The choice of duahty interval 

The sum-rule results for the AB^ depend on the choice of the parameter sq defining the 
upper limit in the duality integrals in (!20|) . (1271) . For the hadronic part the best accuracy 
is obtained by considering small values of sq for which saturation by the ground state is 
a justified approximation. The OPE side, on the other hand, is best suited for inclusive 
quantities for which perturbation theory is valid. The quantity sq must be chosen in 
such a way as to balance between these two cases, and the exact value to use is thus a 
matter of judgement. A useful guide for determining its value is to use QCD sum rules 
for the matrix elements f|T^ to express Sq in terms of the decay constant fs- This makes 
the finite energy sum rule analysis of the three-point correlator parametrically free and 
the analytical results simple, allowing us to discuss qualitative features which are less 
transparent in a purely numerical analysis. 

The two-point sum rule for the decay constant Jb is obtained in the standard way. 
One evaluates the spectral density for the two-point function in both a phenomenological 
hadronic picture and in the OPE. Equating the integrals of the two spectral densities over 
a duality interval gives a result for the decay constant fs- We calculate the OPE spectral 
density by evaluating the two-point function in its crudest approximation, including only 
the bare quark loop. For the two-point function of axial vector currents we have 



Ap^r^iE„E,) 



1 1 r3 

487r2 47r2 [2 



{g^.GG) - 6n^{sGs) (5(^0 + 5(^2)) 



(35) 



z = m\/ s. 



(36) 



For the phenomenological spectral density we have 



p''^{s) = ms-Mi). 



(37) 
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Equating the two expressions as in (|T9l) (local duality finite energy sum rules [19]) yields 

{2nfBY = so(l - zof (AV current) (38) 

zq = ml/ sq. We see that the duality interval parameter sq can be expressed through 
and mf,. We rewrite the expression (l38|) in a form suitable for HQET by substituting 
So = {rrib + Eof and expanding in the ratio Eo/rrib -C 1. Retaining the leading term of 
the expansion we find an equation relating the HQET sum-rule parameter Eq with the 
physical quantity fs- 

{2'KfBf ^ ^El (AV current). (39) 

For = 240 MeV and = 4.8 GeV one finds Eq = 1.1 GeV. Repeating the analysis 
for the pseudoscalar interpolating current, where (neglecting the strange-quark mass) 

(40) 

we have in the HQET limit 

{2TTfBf ^ ^El (PS current) (41) 

which gives Eq = 1.4 GeV. Thus, the numerical value of the duality interval fluctuates 
depending on the channel chosen for its determination. At any rate the results are con- 
sistent with the general expectation that the scale of duality in hadronic physics is about 
1 GeV. 

This idea of determining the value of the duality interval from two-point sum rules 
works well quantitatively also for light quarks. Indeed, by comparison, for light u-, d- 
quarks one finds the relation (27r/^)^ = sq, which gives sq = 0.7 GeV^ for f^, = 130 MeV. 
This is the actual duality parameter for sum rules in the axial-vector channel of light 
mesons PU] . 

The relations (15^ and fHTl) allow for a simple parameter-free analysis in the HQET 
limit. They show the correct scaling for the semileptonic decay constant with the heavy 
quark mass, /b ~ 1/ ^/fn^|, and upon using them in the sum rules the explicit results for the 
bag parameters become independent of rrib, as appropriate for hadronic quantities. For a 
quantitative comparison with full QCD higher-order corrections in EQ/mb are important 
numerically, as the expansion parameter Eo/mb ~ 0.2 is not very small. We see this 
further in our analysis with Borel sum rules. 
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4.2 Finite energy sum rules in HQET 

In this section we present the analysis using the finite energy sum rules ( l20l [2T|) expanded 
to leading order in HQET. We work at leading order in as and ignore even leading-log 
resummation. At this level of precision the HQET approximation can be obtained by first 
evaluating matrix elements in full QCD and then expanding as described in Section 12. 1[ 
To evaluate the phenomeno logical side of the sum rules we use the explicit expressions 
([7]), and to evaluate the OPE side we use the HQET results from the appendix. 

4.2.1 Leading order operators Q and Qs 

We start our analysis with the leading-order operators Q and Qs-, for which we describe 
the procedure in some detail. 

Operator Q with axial vector interpolating current: On the phenomenological 
side of the sum rule ( !20l) we have after subtracting the factorized contribution (cf. eq. [27]) 



where in the second equality we used the HQET limit. To evaluate the OPE side in the 
same limit we use Sj = (m?, + Ei)"^ in the QCD spectral density from the Appendix and 
expand to leading order in Ei/mi,, leaving 

487r^ TT S1S2 

' ^-6). (43) 



487r2 47r2 ml 



Performing the integration on the OPE side we arrive at the sum rule 

|Ai?(27r/B)^ = -2(^,2^2)4- (44) 

Using fl39]) to trade {211/3)'^ for Eq we find the simple result 

AB = _ AM^. (45) 
256 E^ ^ ' 

The result for the non-perturbative bag parameter is independent of m^, as it should be 
in the HQET limit, where dynamical quantities depend on soft physics only. This fact 
can be noticed already from (jUj) by using the scaling relation ~ 1/ ^/Tnb deduced from 
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(|39|). Taking the value of the gluon condensate as {glC^) = 0.48 GeV^ = (0.83 GeV)^ [16J 
we have 

AB = -0.006, (46) 

at ii^o = 1 GeV, which shows that the non-factorizable contribution to the matrix element 
is tiny. 

Examining the expressions for AB, one sees that it is the suppression by the combi- 
nation of variables {27if 3^1711/ Eq = 64_Eq = (3.3 GeV)"^ which leads to this result. This 
combination does not scale with rrib in the HQET limit and Jb is further enhanced by 
in the laige-Nc limit. Since the scale of the gluon condensate is given by (0.83 GeV)*^, 
the result for AB is proportional to the fourth power of a small number. In the absence 
of any accidental numerical enhancement of the coefficients, which we do not see, the 
"natural" size of the deviations from factorization is extremely small. 

The answer (H6l) is the leading-order HQET result. To get a feel for the size of the 
subleading terms, we list the next few terms in the expansion of the OPE spectral density: 



AB 



3 




256 




3 




256 




3 




256 


El 



3 mi) KrUhJ 9 KniiyJ 
[1-0.8 + 0.4-0.2 + ...] 

(0.5) , (47) 

where we used mb = 4.8 GeV and to obtain the last line we evaluated the full QCD 
result. This shows that the subleading terms are not small, and that keeping only the 
leading-order term misses the full QCD result by a factor of two, at least at i?o = 1 GeV. 
Given the small size of AB the factor of two is numerically irrelevant, and is actually 
within the uncertainties of the analysis. We return to this point in Section El using the 
subleading operator R2 as an additional example. 

Operator Q with pseudoscalar interpolating current: We can repeat the compu- 
tation using a pseudoscalar interpolating current. At leading order in l/m^ we neglect 
rris and expand as before, finding 

^AE(27r/B)^ = i-{glG')El + M'{-sGs)E,) . (48) 

The mixed quark-gluon condensate is parameterized as (sGs) = ml{ss). For numerical 
evaluation we use = 0.8 GeV^ [211 [22] and (ss) = 0.8 (mm) [231 [2l]. For the light quark 
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condensate {uu) we take {uu) = (—0.24 GeV)'^. A convenient normalization for the mixed 
quark-gluon condensate is 

TT^sGs) = -0.1 GeV^ = (-0.63 GeV)^ (49) 

which is of the same order of magnitude as {glG"^) and is really given by the hadronic 
scale of 1 Ge V (or say by the p- meson mass mp = 770 Me V) . We see that all dimensionful 
quantities are on the order of the fundamental QCD scale of 1 GeV as expected. 
Using (!39|) to eliminate fs one finds 

3 fig^^G^) A7T^(sGs)\ 3 , 



at i?o = 1 GeV. The result is more or less the same as with the axial-vector interpolating 
current but the structure of contributions changed. The gluon condensate contributes less 
and mixed quark-gluon gives a contribution (it was zero for the axial- vector case). 
Using Eq^ = 1.4 GeV from fHT!) for the pseudoscalar channel one finds 

^B = --( - ^^^^^] = -- (0.062 + 0.074) = -0.006 (51) 
64 V2(^o^'^)^ (^o^'^)^ ; 64 ^ ' ^ ' 

which coincides with the previous result from the AV current. Even though it is more 

proper to use the result fj4T|) for the PS channel, for the remaining operators we shall use 

fl39|) in both the AV and PS channel; the differences are very small. 



Operator Qs with axial vector interpolating current: Repeating the analysis for 
Qs with an axial-vector current we find the sum rule 

- ^ABs{27rfsy = ^ [l{9'sG')E', - A7c\sGs)Eo) (52) 

and 

ABs = -— ( - 247i'^-^] = -— (0.48 + 2.4) = -0.005 (53) 

* 640 V E^ El y 640 ^ ' ^ ' 

again a very small number. Note that the result is dominated by the contribution of the 
mixed quark-gluon condensate, even though it is formally subleading compared to the 
gluon condensate. This is a general situation, as contributions from the gluon condensate 
are often small numerically. 



17 



Operator Qs with pseudoscalar interpolating current: Repeating for the pseu- 
doscalar current we have 

640 V ^0 ^ ^0 / 640 



^Bs = I - Svr^l-^ 1 = -— (0.48 + 0.8) = -0.006. (54) 



The coefficients of the gluon and mixed quark-gluon condensates changed, but their sum 
is very close to that obtained with the axial-vector current. 

We conclude that the deviation from factorization is tiny just because of the scales 
involved. No surprisingly big numbers or drastic cancellations occurred in the analysis. 

4.2.2 Subleading operators R2 and 

The analysis is essentially unchanged for the subleading operators R2 and R^. The new 
feature is the appearance of the parameter A = Mb — rn^ even at leading-order in the 
HQET expansion (numerically Mb = 5367.5 ± 1.8 MeV [25j). It enters through the 
expansion of the matrix elements ([7j), which read 

(5|i?2|5)HQET = -flml [2Km,) [l - ^) B2 (55) 

(S|i?3|5)HQET = flml[2Km,)[l + ^^B,. (56) 

This A/mf, power suppression of the matrix elements on the phenomenological side of 
the sum rules is compensated by an Eo/rrib suppression from the OPE spectral densities. 
Then, up to a factor of i?o/2A ~ 1, the magnitude of AB for the subleading operators is 
fixed by the same scale ratios as before, and as with the leading-order case there are no 
large deviations from factorization. 



2Am,(27r/B)'^Ai?2 = ^ (^{g'',G^)m,E', (^-^) - 2n^{rsGs)m,E^'^ (57) 



Operator R2 with axial vector interpolating current: 

I 2 . „ 1 / , n — ^-Ov / 3 

Notice that for this subleading operator the phenomenological and OPE sides of the sum 
rule are suppressed by the hadronic scales A and Eq respectively. This is explicit in the 
HQET expressions but not in the QCD ones. Using fl39l) and taking A = 0.5 GeV 

Aft = ^ A f + lr(£££>) . A (0.48 - 0.13) = 0.012. (58) 

2A 256 V E^ 3EI J 256 ^ ' ^ ' 

We see that AB2 is again very small, although this time it is a positive number instead 
of a negative one. 
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Operator R2 with pseudoscalar interpolating current: 



Operator R3 with the axial vector interpolating current: 



(59) 



and 



7, Eo 1 ( (glC^) fl\ 2tt\sGs)\ 1 , 



A53 = 0.004 (61) 



Operator with pseudoscalar interpolating current: 



-A53 = ^— ' = — 0.12 + 0.1 = 0.003 62 



and 



A53 = 0.003. (63) 

We can summarize by saying that the finite energy sum rules within the HQET ap- 
proximation suggest that factorization is perfectly precise, if only non-perturbative con- 
densate effects are taken into account. The bag parameter Bi has the largest violation of 
factorization, but it is still very small in absolute terms, approximately 1%. 



5 The bag parameters from Borel sum rules 

We have seen in the previous section that the deviations from factorization for both 
the leading and subleading operators are very small. In this section we perform a more 
thorough numerical analysis using Borel sum rules. This serves to confirm these results 
and to show that it is possible to impose very conservative error estimates without altering 
this conclusion. We also use the numerics to compare the HQET and full QCD results. 

5.1 Borel sum rules in full QCD 

The Borel sum rules in full QCD are evaluated according to (|22l) and the analogous 
expression for the pseudoscalar interpolating current. Although it is possible to evaluate 
the Borel integrals analytically, the results are quite lengthy and we do not need them for 
this purely numerical analysis. 
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Figure 3: Plot of —AB vs. (GeV^) for the leading-order operators Q and Qs and AB 
for the subleading operators R2 and -R3 with the Borel sum rules in full QCD. The short- 
dashed lines are obtained using an axial vector interpolating current, and the long-dashed 
lines using a pseudoscalar current. The parameter values are given by = 4.2 GeV, 
fB, = 240 MeV, So = 36 GeV^ and TUg as explained in the text. 

To evaluate the sum rules, we must first give numerical values for the QCD parameters 
/ss) ^b, and rris. For the decay constant we choose Jb^ = 240 MeV as the default value. 
For the 6-quark mass one can take the pole mass or the MS mass. The pole mass is 
mf^" = 4.8 GeV while the MS value is mf^ = 4.2 GeV [26l [27|. For the full QCD 
analysis the MS mass is more appropriate. However, since we are working to lowest order 
in as, we cannot distinguish these two quark- mass definitions, and the difference can be 
accounted for as an additional uncertainty in AB. This difference would be under control 
if as corrections were taken into account. 

The strange-quark mass appears on the OPE side of the sum rules for all channels, 
and in the phenomenological side for the case of the pseudoscalar interpolating current. 
We have seen that AB is extremely small in all cases, and the effects of a non-zero strange 
quark mass do little to alter this. We choose to keep it non-zero on the phenomenological 
side for the leading-order operators Q and Qs, using nis = 100 MeV. Keeping it non-zero 
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Figure 4: Plots of {—)AB vs. W (GeV) obtained with the Borel sum rules in HQET. 
The short-dashed lines are obtained using an axial vector interpolating current, and the 
long-dashed lines using a pseudoscalar current. We take A = 600 MeV, nib = 4.8 GeV, 
fB, = 240 MeV, and 2Eo = 2.5 GeV. 



in the OPE spectral densities complicates the analytical expressions without changing the 
final results in a significant way. For the subleading operators it is consistent to set it to 
zero at this order in the heavy-quark expansion. 

The results for AB vs. the Borel parameter for each operator are shown in 
Figure [31 The two lines in each plot are obtained by using an axial-vector and pseu- 
doscalar interpolating current. All results have a reasonable stability region in at 
10 GeV^ < < 20 GeVl There is some dependence on the choice of interpolating 
current, which as we will see later is within the uncertainties of the analysis. We note 
again that AB is positive for the subleading operators R2 and -R3, whereas it is negative 
for the leading operators Q and Qs- 



5.2 Borel sum rules in HQET 

The Borel sum rules in HQET are performed according to (!28|) and analogously for the 
PS interpolating current. We focus on numerical results, although the analytical results 
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Operator 


AB(%) QCD 


AB(%) HQET 


Q 


-0.6 ±0.5 


-0.6 ±0.5 


Qs 


-0.5 ±0.4 


-0.6 ±0.4 


R2 


0.3 ±0.3 


0.8 ±0.7 


R3 


0.3 ±0.2 


0.3 ±0.2 



Table 1: A summary of the results. 



for the Borel integrals are very simple. In fact, they reduce to those from the finite energy 
sum rules in the limit W ^ 00. In contrast to our treatment of the finite energy sum 
rules, however, in our numerical studies we treat Eq and f^^b as free parameters. We 
again use fs^ = 240 MeV as the default value. While in the QCD calculation the MS mass 
was more natural, in HQET the pole mass appears in the construction of the effective 
theory and is more natural. We use m^"'*^ = 4.8 GeV. 

The results for AB vs. the HQET Borel parameter W for each operator are shown 
in Figure m The plots are stable in the region 1 GeV < W < 2.5 GeV, which is rather 
typical for Borel sum rules in HQET. The values of AB in the stability range are close 
to those in the QCD plots in Figure [3l The one noticeable exception is R2, where the 
HQET values are about twice as large as the QCD ones. We comment further on this in 
the next section. 

6 Final results and discussion 

We now present our final numerical results and estimate the associated uncertainties. The 
results are summarized in Table [H 

To obtain the table entries for full QCD, we fix the Borel parameter at = 15 GeV^ 
and vary the other parameters in the ranges 

QCD 

210 MeV < /b, < 270 MeV 

4 GeV <mb< 4.4 GeV 

32 GeV^ < So < 40 GeV^ 

where the default values lie in the center of the above ranges. We also vary the condensates 
about their default values by ±30%. For a given case, we find upper and lower values 
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-AB(%) /b variation (HQET) -AB (%) A variation (HQET) 




Figure 5: Plots showing — AB vs. W (GeV) for the case of the Qs operator, axial-vector 
current, using a range of parameter values in the HQET sum rule. The plots are made 
by varying the parameters as explained in the text. For the condensate variations, the 
dark-gray band corresponds to the gluon condensate and the larger light-gray band to 
the quark-gluon condensate variation. In each case, the dashed line corresponds to the 
lower border of the parameter range (e.g. /b^ — 210 MeV) and the solid hue to the upper 
border (e.g. fe, = 270 MeV). 

of AB to identify the error ranges. For the /b^ , and Sq variations the ranges are 
asymmetric; in those cases we use the larger deviation in the error analysis. Finally, 
we add the uncertainties from each of the five variations in quadrature, and average the 
results from the axial- vector and pseudoscalar interpolating currents to obtain the results 
quoted in the table. 

The procedure is the same for the HQET sum rules, although the set of parameters 
is different. This time we fix the Borel parameter ai W — 2 GeV and the 6-quark mass 
at rrib — 4.8 GeV, and vary the other parameters in the ranges 

HQET 

210 MeV < fe, < 270 MeV 
500 MeV < A < 700 MeV 
1 GeV <Eo< 1.5 GeV 
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where the default values lie in the center of the above ranges. The condensates are again 
varied by ±30% about their default values. The final table entries are obtained as for the 
QCD case. 

To illustrate the uncertainty associated with each parameter variation, we choose as 
an example the Qs operator with an axial- vector interpolating current in HQET. The 
range of AB associated with each variation is represented by the gray bands in Figure [51 
It is seen that the largest errors are associated with the value of the decay constant fs- 
This is not surprising, since the explicit results scale as 1//^. At the default value W = 2 
the dependence on A and Eq is moderate. The results depend linearly on the condensates 
and at W = 2 the uncertainty due to the condensates is comparable with that due to /b- 

In all cases except for R2, our central values for AB in QCD and HQET turned out 
to be (nearly) equal. However, in interpreting this result, one should be clear that not 
only the bag parameters, but also the QCD parameters Jb and Mb have an expansion 
in l/ruf,. When comparing the QCD result with the HQET result, we have no means of 
disentangling the corrections to fs and Mb from those to Bi, so it is not obvious whether 
numerical discrepancies are due to corrections to the bag parameters, form factors, meson 
masses, the OPE, or even our choices of sum rule parameters. The conclusion to make is 
that the leading-order expansion and the full results are consistent with one another in 
all cases, within the uncertainties of the analysis. 

This said, further investigation of the HQET series for the OPE spectral densities for 
R2 reveals some interesting features. As an example, we take the piece of the spectral 
density for R2 multiplying (GG) as calculated with an axial-vector current, and consider 
some higher-order terms in the Eq/itli, expansion of the integrated spectral density. Using 
the notation Xi = Ei/rrib, integrating over the square < Xi < Xq, and normalizing to the 
leading-order term in the xq = Eq/uii, expansion, we have 

4/dx,dx.p-(x„x.) 154^„ + ^,^ + ... = 1.0-1.2 + 0.9 + ... (64) 



4/da;ic/x2(-9/2xi -9/2x2) 27 27 

To derive the numbers we used xq = Eo/rrib ~ 0.2 for Eq = 1 GeV. The second and third 
terms are as large as the first, and the corrections do not fall below 10% until the sixth 
term, so the "HQET" expansion is not well behaved. We put HQET in quotes, because 
the expansion is just the diagrammatic one, not a rigorous one in terms of operators. It 
would be interesting to see whether this poor convergence persists even with a more careful 
treatment of the subleading corrections. If so, this would have important implications for 
lattice QCD results, where corrections to the HQET limit are not easy to control. 

In quoting our final results, we used only those obtained from the Borel sum rules. 
However, one can work with either finite energy or Borel sum rules. Finite energy sum 
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rules can be obtained from Borel sum rules in the limit —>■ oo and are therefore more 
sensitive to the model of the continuum. We used both and saw little difference. Our sum 
rule analysis is by no means unique. For instance, one can change the duality integrals by 
modifying each side of the sum rule in the same way (for instance by dividing both sides 
by (S1S2)). This definitely changes the shape of the curves and can provide better stability. 
However, our main point is that AB is so small that we need not be too sophisticated with 
the sum rules analysis. The splitting into factorized and non-factorized parts is powerful 
and useful precisely because the absolute value of AB turns out to be small. Even with 
very conservative error estimates the results are numerically informative, and our final 
results - the range for the values of AB - rather reliable. 

It is instructive to compare our approach to lattice QCD. In the lattice approach 
the parameter B is computed as a whole, since a splitting into factorizable and non- 
factorizable parts is not possible at the level of simulation. Then for the computation of the 
parameter B (and not AB directly) even good accuracy of the method (say, about 20%, 
a typical accuracy in hadronic physics) gives a less precise statement about factorization 
than our technique. 

Our analysis was limited to leading order in perturbative corrections. A more accurate 
determination would require the computation of the next-to-leading order perturbative 
contributions. These involve three-loop diagrams and this is a non-trivial task. Results 
are nonetheless available for the operator Q [18], where it was shown that these corrections 
amount to about 10%. For the other operators, we can say only that the corrections are 
parametrically on the order of ag/vr and are also expected to be around 10%. Thus, a 
qualitative prediction of the sum-rule analysis is that deviations from factorization are 
suppressed either by scale ratios or by the strong-coupling constant and are therefore 
small. 
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7 Conclusions 



We used QCD sum rules to calculate the bag parameters for the leading and next-to- 
leading order operators in the l/m;, expansion of the transition operator used to analyze 
Bg — Bg mixing. We found that the violation of the factorization approximation for the 
matrix elements of both the leading and subleading operators due to non-perturbative 
vacuum condensate contributions is well under control and small. Our final results for 
the parameters Ai^j = Bi — 1 are 

AS|q = -0.6 ±0.5% 

AS|Qg = -0.5 ±0.4% 

A5|^2 = 0.3 ±0.3% 

AB|fl3 = 0.3 ±0.2% 

We believe that our very conservative error estimates make our final quantitative results 
(the range for the values of A5) rather rehable. 

Our result that the non-perturbative contributions to the ASj are extremely small 
means that the non-factorizable contributions to the matrix elements are most likely 
dominated by calculable perturbative effects. We did not attempt to include the next-to- 
leading order perturbative contributions in our analysis. Our naive expectation, based on 
the existing calculations for the operator Q, is that these corrections can contribute an 
additional ±10%. To clarify this point would require the evaluation of the set of three- loop 
diagrams appearing in the perturbative analysis. 
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SI 349/10-1). BDP acknowledges the support of the SFB/TR09 "Computational Particle 
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8 Appendix 

Here we compile the results for the condensate contributions to the OPE spectral density 
in each of the eight cases. For the axial- vector current we single out the scalar amplitude 
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multiplying the structure tensor structure P1P2 ■ For the pseudoscalar interpolating current 
there is only one amplitude as the correlation function is a scalar. 



8.1 Spectral densities for Q 

For the AV interpolating current we have 

ApAv -- 



I as 1 

{—GG) {piP2)2ZiZ2{-3 + Zi + Z2- 2Z1Z2) 



487r^ TT ■ S1S2 

1 1 1 



487r2 47r2 ml 



(65) 
(66) 



where z — ml/s and we omit the factor 9{si — ml) setting the lower limits of integration 
for the Si. The factor (P1P2) — si/2 + S2/2 ior — 0. To take the heavy-quark hmit in 
the second line we used Si — m^(l + Xi)^ and expanded to lowest order in Xi — Ei/mj,. 
For the PS interpolating current: 

1 X 

^PPS = 7^(— G'C) {piP2)ml3{-2 + z, + Z2-ZiZ2) (67) 

+ {PiP2)-r^{sGs)mb (—2ziS{s2 - ml) + —2z25{si - ml)) (68) 

- 2>{glG'^) + 12m67r^(sGs) (25{s2 - ml) + 25{si - ml)) (69) 



487r2 47r2 



8.2 Spectral densities for Qs 

AV interpolating current: 
1 



GG) 



ml 1 



(6 - 3(2:1 + Z2) + Z1Z2) 



487r2 ^ -K ' S1S2 2 
+ T^(^<^«)^& f-(-2 + Zi)Ks2 - ml) + -(-2 + Z2)5{s^ - ml) 

-{glGG) - l2m,Ti\sGs) (5{s2 - ml) + 5{s^ - ml 



1 1 1 
487r2 47r2 ^ 



PS interpolating current: 
Apps = 



(6 - 82:1 - 82:2 + 2:12:2) + {PlP2f'ZiZ^ 

^{sGs)mb ((-2 + ^i)(5(s2 - ml) + (-2 + Z2)5{si - ml)) 



-^l^GG)^ 



167r2 
1 1 



487r2 47r2 



■{glGG) - 12mb7i^{sGs) (6{s2 - ml) + S{si - m, 



(70) 
(71) 



(72) 
(73) 
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8.3 Spectral densities for R2 

AV interpolating current: 



1 ,a 



AV 



487r2 TT 

9ziZ2 — 2zi + 1 



— {-^zlzl + I2zlzl - Azlz2 - 9zlzl - ?>z\z2 



■z\z2{2Z]_z\ — 4:ZiZ2 + Zi+ AZ2 — 3) 



+ {Zl Z2) 



+ 



S1S2 



16/1^ 

1 1 1 



487r2 47r2 



9 



xi{glGG) — 12ximbn'^ {sGs)S{s2 — ml) + {xi <-> X2) 



(74) 
(75) 



Note that the l/m^ suppression of R2 compared to Q and Qs is manifest only after the 
HQET expansion. The coefficient of the (g^GG) term is large and there is a relative sign 
of {GG) and {sGs) terms compared to all other cases. This is a unique feature. 
PS interpolating current: 

■3 



ApPs = T^(-GG) 
487r^ 71 



_^3 (piP2) ^^^^^^2 _ 32;^^2 + -Zi + 3^2 - 2) 



2 S1S2 
1 



— ^1^2 ~^ 3z^Z2 — -2^ — 22^1 2^2 — -^1 ~l~ 2) 
+ (^1 ^ 2^2) 



iDTT^ 777.6 
1 1 



'!lii4-l) + (l^,li-,, + l)]6is2-rnl] 



487r2 47r2 



xi{ggGG) + 12a;imb7r {sGs)6{s2 — m^) + (xi 0:2) 



(77) 



8.4 Spectral densities for i?3 

AV interpolating current: 



AV 



^ - 4siS2 



487r2 ' TT 



2:1 (-2:1 2:2 + 72:1 2:2 - 32;i - 92:2 + 6) 



+ {zi ^ Z2) 



:{sGs) 



{P1P2) ^ 



-z\ + 3^1 - 2)5(s2 - fn^) + {zi ^ Z2) 



111 

487r2 47r2 ^ 



-xi{glGG) — 12ximb7r'^ {sGs)S{s2 — ml) + {xi ^ X2) 



(78) 
(79) 



PS interpolating current: 

487r^ TT 



(P1P2 ) 
16 



{-Ztz2 + 10ziz2 - 6Z1Z2 - 12zi + 9) 



28 



(PlP2)^ 2 / , in1 , / ^ 

■— ZiZ2{-Z2 + 1) + [Zl ^ Z2) 

-^3182 



+ ^{sGs)^-^(-zf + 3z, - 2)S{S2 - ml) + (z, ^ Z2) (80) 
1 1 



3 

-^Xi{glGG) — 6ximh7i'^{sGs)6{s2 — ml) 



+ [Xi ^ X2) 81 



487r2 47r2 

8.5 Four-quark condensates in the OPE 

In this subsection we discuss the treatment of the 4-quark condensates, and show that their 
non-factorizable contributions to the current correlator vanish upon applying the vacuum 
saturation approximation. The Fourier-transformed three-point correlators involving Q 
and Qs are given by 

Tip,,p2) = 2{0\{sTj{mh+]^i)Tos){sTjim,-]^2)ros)\0) (—^ ^] , (82) 

\ ~ Pijyi^b ~ P2) / 

where Tq = 7^(1 — 75) for Q and (1 — 75) for Qs, and Tj = 7*^75 for the axial-vector 
interpolating current, and 275 for the pseudoscalar one. In contrast to the case for the 
gluon and quark gluon-condensates, explicit expressions for the factorizable contributions 
to the three-point functions are needed. They read 

Tf,, = 2 (1 + i-) n(p0n(-P2) (83) 

for Q and 

r?ac = 2 f 1 - TT^) Il{p^)U{-p2) (84) 



for Qs- The functions Il{pi) for an operator with Dirac structure Tq and an interpolating 
current structure T j read 

U{p) = I j rfW^'^(0|T{srj6(x)6ros(0)} |0) = {{)\sTj{m, + ^)Tos\0) . / ^. . (85) 
J ^ J [mi — p^) 

We now isolate the part of the current correlator AT which contributes to a non-zero AB 
by using the definition (|23|) . which gives 

Ar(pi,P2) = T\p,,p2) - 2 (1 + ^) n(pi)n(-P2). (86) 

The full correlator T contains four-quark matrix elements of the form (srissr2s). Eval- 
uating these matrix elements using (132|) . we find that AT vanishes for all cases. 
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